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Abstract 

Sixteen and thirty-two ply quasi-isotropic laminates fabricated from 
AS4/3501-6 were subjected to pure tension, pure torsion, simultaneous 
tension and torsion, and torsion fatigue. Layups tested were 
[45n/-45 n /0n/90 n ]s, with n=2 or 4. A torsion damage pattern consisting 
of a localized matrix crack and delaminations was characterized, and the 
measured torsional stiffnesses were compared with calculated values. It 
was found that a combination of tension and torsion led to failure at 
smaller loads than either type of deformation acting alone. Further work 
is required to determine the exact form of the failure criterion. 

Introduction 

The delamination of multidirectional composite laminates subjected 
to tensile loads is fairly well understood. The effects of lay-up, moisture, 
thermal residual stresses, and resin toughness have been described. 1 ’ 2 In 
some structural applications, including helicopter rotors, composites are 
subjected to a combination of tension and torsion. In this report the 
damage induced in quasi-isotropic graphite/epoxy laminates subjected to 
pure torsion and combined tension and torsion loads will be described and 
contrasted to the damage observed in these same laminates when 
subjected to pure tension loading. 

Experimental Procedures 

Sixteen and thirty-two ply quasi-isotropic laminate plates with a 
[45 n /-45 n /0n/90 n ]s layup (where n= 2 and 4) were fabricated from 
Hercules AS4/3501-6 graphite/epoxy prepreg using the manufacturer’s 
recommended autoclave cycle. Specimens 1" wide and 5" or 10" long were 
cut from the plates. 

Tension-only tests were run in stroke control, using a clip gauge to 
monitor strain. Tension/torsion testing was performed using a 
servohydraulic test frame equipped with hydraulic wedge grips. Torsion 
tests and combined tension/torsion tests were run in rotary stroke control 
and axial load control. Angles of twist were determined by the rotary 
transducer mounted on the hydraulic actuator. Axial load was monitored 
and verified to be constant and independent of the angle of twist. Static 
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tests were run at a rate of 0.125 degrees/sec; fatigue experiments were 
performed with a load ratio of R— 0.10 under sine-wave loading at 5 Hz. 
Specimen edges were coated with water-based typewriter correction fluid 
to aid in visual detection of the onset of damage. Critical strains and 
critical angles of twist were recorded when delamination occured. After 
testing, some specimens were treated with a dye-penetrant that is opaque 
to x-rays, then x-rayed. Specimen edges were later polished for 
photomicroscopy. 

F.xperimental Results 


Critical delamination strains, critical loads, and moduli obtained from 
static tests in pure tension are summarized in Table I. Failures in this type 
of tensile specimen have been described previously. ^ Delaminations form 
in the 0/90 interface, jumping through the 90° plies to the symmetric 90/0 
interface randomly along the entire length of the free edge between the 
grips. 

The data for pure torsion are shown in Table II. The critical 
normalized angle of twist for delamination, 9/L, is relatively independent 
of the gauge length but decreases with increasing specimen thickness. 

The damage pattern observed on the edges of the 16 and 32 ply 
specimens is shown in Fig. 1. On each edge, one oblique crack in the 
central 90° plies has a pair of delaminations growing at either end of it in 
opposite directions. This "primary crack is presumed to be the point 
where the delaminations initiated. A schematic of the characteristic failure 
mode is shown in Fig. 2. The delaminated areas show up quite clearly in 
the dye penetrant enhanced radiographs. The hourglass-shaped shadow in 
Fig. 3 marks the boundaries of the edge delaminations. Furthermore, both 
the x-ray and optical microscopy (Fig. 4) reveal the primary crack as well 
as other secondary matrix cracks that develop in the 90° plies. 

The effects of combined tension/torsion may be studied using the 
same type of specimen. To do this, a constant axial load, equal to 
approximately 1/2 the failure load in pure tension, was applied to the 
specimen as the torsion load was introduced. Figure 5 shows the 
delamination onset data from the pure tension, pure torsion, and combined 
tension and torsion tests, with the axial strain shown on the ordinate and 
the normalized angle of twist shown on the abcissa. Curves are faired 
through the maximum and minimum values to bound the data. Figure 6 
shows the delamination data in terms of the axial load and torque. 

Based on the elasticity solution for the torsion of a rectangular 
prismatic bar, 4 the maximum shear stress should occur at the midplane of 
the laminate thickness. Because of the laminated nature of composites, 
however, deviations from pure shear may occur near the edge. The matrix 
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crack visible in Fig. 4a. makes an angle of approximately 50° with the 
specimen face. Among all the specimens tested in pure torsion this angle 
ranged from 37° to 65°, with a mean of 51°. The corresponding cracks in a 
tension specimen lie at 90° (perpendicular to the tensile direction.) It was 
therefore anticipated that as greater proportions of tension loading were 
added in tension-torsion tests, the angle would increase toward 90°. The 
data from all the micrographs examined are collected in Table III. As had 
been expected, delaminations did form at a smaller angle of twist in these 
combined loading specimens than they had in specimens tested in pure 
torsion. However, the matrix cracks made an average angle of 52°, only a 
slight deviation from the pure torsion result. 

Preliminary torsion fatigue data are displayed in Fig. 7. The logarithm 
of the number of cycles to delamination onset increases steadily with 
decreasing strain. No critical strain threshold (leveling out of the curve) is 
apparent out to 10 6 cycles. 

Estimation of Torsional Rigidity 

From strength of materials 6 , the torsional rigidity of an isotropic 
prismatic bar is defined as the product of the shear modulus, G, and the 
polar moment of inertia, J, and is determined as 

JC=IL 

e (1) 

where T is the applied torque, L is the length of the bar, and 0 is the 

angular deformation of the bar measured in radians. In plate theory 6 , the 
torque is given by 


T = M xy w (2) 

where M xy is the twisting moment per unit width, w, and the twisting 
curvature is given by 


ic -A. 

K xy - L 

Substituting Eq.(2) and Eq.(3) into Eq.(l) yields 

JG = ^L= MxyW 
0 K *y 


(3) 


(4) 
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In laminated plate theory?, the twisting curvature may be related to the 
twisting moment by the following 

Kxy = D66 ^ M xy (^) 


where D66 1 represents the last term in the inverse bending stiffness 
matrix. Substituting Eq.(5) into Eq.(4) yields an estimate for the torsional 

rigidity of a composite plate as 


JG = 


TL _ w 

e D'i 


( 6 ) 


Therefore, for a composite laminate of known width, lamina properties, 
stacking sequence, and ply thicknesses, we can use Eq.(6) to estimate t e 
torsional rigidity of the plate. Equation (5) assumes that a torque is 
applied to all four sides of the laminated plate. In our study, however, the 
two longitudinal sides are stress free. Hence, laminated plate theory will 
underestimate the torsional rigidity because of the discrepancy between 
the actual boundary conditions and those assumed implicitly m Eq.(b). 

A more accurate estimate of the torsional rigidty for these quasi- 
isotropic plates may be obtained from the elasticity solution for an 
isotropic prismatic bar of rectangular cross section loaded in torsion at the 
ends. 4 The torsional rigidity from this solution may be expressed as 


JG = — = wt 3 G(J 
0 


(7) 


where 


P = Mi -»•«£) 


and t is the plate thickness and G is the shear modulus. This solution now 
incorporates the appropriate boundary conditions. However, the laminate 
is assumed to be isotropic. The quasiisotropic composite laminate is 
isotropic in the plane of the laminate, but the through-thickness properties 
are anisotropic. Hence, G xy = G xz , but G xy is not equal to G yz . T ere ore, e 
elasticity solution for an orthotropic prismatic bar of rectangular cross 
section loaded in torsion at the ends? may be more appropriate. The 
torsional rigidity from this solution may be expressed as 
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jG=IL =wt 3 Gx/ p(c) 

0 ( 8 ) 


where 



and 



The properties in table V were assumed in order to calculate G xy 
from laminated plate theory as 


Gxy — 



(9) 


where A 66 is the last diagonal term in the inverse extensional stiffness 
matrix. To determine the torsional rigidity, G was set equal to G xy in Eq.(7) 
and G xy = G xz in Eq.(8). Two different values of G yz , the out-of-plane shear 
modulus, were reported for AS4/3501 -6. in Reference. 8. The value of G yz 
measured for an orthotropic [0/45/-45/0h s laminate was less than G yz 
measured for a unidirectional laminate, illustrating the anisotropic nature 
of this property. No value was available in the literature for a quasi- 
isotropic laminate similar to the ones tested in the present study, so both 
values of G yz from Table 2 of Ref.[8] were used in Eq.(8) to predict the 
torsional rigidity of a quasi-isotropic laminate. Only the first 5 terms (n = 
1,3, 5, 7, 9) in the series expression for p in Eq.(8) were required for this 
series to reach a converged value. 

Figure 8 shows the experimental values of torsional rigidity 
measured in this study (EXP) as well as calculated values. The bars labeled 
LPT were estimated from laminated plate theory using Eq.(6). The 
isotropic elasticity results from Eq.7 are designated ISO and the two sets of 
orthotropic elasticity results from Eq. 8 are designated ORTHO ELAS. 
Comparing these results, we see that the torsional rigidity estimated from 
laminated plate theory using Eq.(6) is significantly lower than the 
experimental values due to the error in implicit boundary conditions noted 
previously. The measured torsional rigidity may not agree precisely with 
Eqs.(7) and (8) because experimentally the angle of twist, 0, was measured 
using the angular deflection of the lower grip relative to the upper grip. 

The assumption of uniform warping used to develop these solutions is 
violated unless 0 is measured between two points on the specimen away 


from the grips. A better measurement of 0 would be expected to lead to 
better agreement with the calculation. Nevertheless, the agreement 
between the experimentally measured torsional rigidity and the 
calculations from Eqs.(7) and (8) is reasonable. 


Discussion 

The effect of combining tension and torsion loads is seen in Fig. 5 and 
Fig. 6.. The combination of torsion and tension results in an earlier 
delamination failure than would be observed by either loading acting 
alone. The data generated in this study indicate that a continuous failure 
criterion may be developed to account for the combination of loads. It is 
not yet clear what parameter should be used for this criterion or what 
form the criterion would have. It is desirable to cast the criterion in terms 
of a generic parameter, one that is independent of structural parameters 
such as layup, stacking sequence, thickness, gage length, etc. The data in 
Figures 5 and 6 clearly show a thickness dependence. Hence, critical 
values of loads (torques) and strains (twist angle per unit length) are 
structural properties, not generic material properties. 

For tension loading, the delamination failure criterion was cast in 
terms of a critical strain energy release rate, Gc, assuming an edge 
delamination geometry!. For torsion or combined tension/torsion loading, 
however, a local delamination model may be needed, similar to a G solution 
derived previously for tension^ because the delamination appears to 
emanate from a dominant matrix crack and grow in an antisymmetric 
pattern (Fig.2). The torsional rigidity terms in a local delamination "G" 
solution for torsion cannot be represented using laminated plate theory, as 
was done in tension, however, because of the discrepancy in edge 
boundary conditions between the test specimen loaded at the two ends 
only and the laminated plate theory assumption of torsion applied to all 
four sides. Further work is needed to resolve these difficulties. 

Other investigators are developing 3D finite element analyses for the 
interlaminar stresses and strain energy release rates associated with local 
delamination induced by torsion or combined tension/torsion loading! 1. 
These models should help explain the local stress states through the 90 
degree ply thickness that result in the observed matrix crack angles (Table 
III). Also, these models will provide benchmark solutions with which to 
calibrate any simple closed-form models that may be developed for strain 
energy release rates associated with local delamination induced by torsion 
or combined tension/torsion loading. 


Conclusions 


Based on the experiments conducted in this study, the following 
conclusions have been reached: 

1. For torsion or combined tension/torsion loading, 
delamination appears to emanate from a dominant matrix 
crack and grow in an antisymmetric pattern. This pattern is 
markedly different from the edge delamination pattern 
observed for tension loading on the same laminate. 

2. For pure torsion loading, the angle of the dominant 
matrix crack with the laminate axis is not simply 45 degrees as 
anticipated based on the maximum shear obtained from the 
elasticity solution for a prismatic bar with a rectangular cross 
section. Furthermore, this crack angle changes very little when 
tension loading is combined with torsion loading, evidently due 
to the influence of anisotropy and the free edge. 

3. Critical values of loads (torques) and strains (twist angle 
per unit length) depend on specimen thickness, and hence, are 
structural properties, not generic material properties. 

4. The addition of torsion to tension results in an earlier 
delamination failure than would be observed under either 
loading acting alone. 

5. Torsional rigidity cannot be represented using laminated 
plate theory because of the discrepancy in boundary conditions 
between the test geometry, which involves loading the 
specimen at the two ends only, and the laminated plate theory 
assumption of torsion applied to all four sides. 
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Tension Edge Delamination Results 
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Table III 

Torsion Crack Angles 
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Gauge 

Length, Angle, degrees 



Table IV 

Tension/Torsion Crack Angles 


Gauge Angle, degrees 



TABLE V 


AS4/3501-6 GRAPHITE EPOXY 
LAMINA PROPERTIES 


En = 19.5 x 10 6 Psi 
E 22 = 1.48 xIO 6 Psi 

G 12 =0.8x10 6 Psi 

V 12 =0.3 


h = 0.005 in. (ply thickness) 




Fig. 1. Photograph of edge of failed specimen showing pattern of 
cracking in brittle white coating 



Fig. 2. Schematic of torsion failure mode 
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[45 2 /-452/02/90 2 3s SPECIMEN FAILED IN 



Fig. 3. Dye/penetrant enhanced x-ray radiograph of failed specimen, 
showing delaminations and ply cracks 



(a) Primary matrix crack and (b) Secondary matrix cracks 

local delamination 


Fig. 4. Optical micrographs of polished edge of 16-ply specimen 
failed in pure torsion 

a) Primary crack, with delaminations growing from both 

ends in opposite directions 

b) Secondary cracks and adjacent single delamination 
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Fig. 6. 



Combined tension/torsion data plotted as critical axial strain 
vs. critical angle of twist 



Combined tension/torsion data plotted as critical axial load vs 

critical torque 
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Fig. 7. Critical Torsion angle as a function of number of fatigue cycles 
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Fig. 8. Measured and calculated torsional rigidities (as described in text) 
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